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ABSTRACT 



The recent report of X-ray Fe features in the afterglow of the gamma-ray burst 
GRB 991216 may provide important clues for identifying the nature of its progenitor 
and constraining the burst mechanism. We argue that the strong line emission can 
be attributed to the interaction of a continuing (but decaying) post-burst relativistic 



CNJ ■ outflow from the central engine with the progenitor stellar envelope at distances less 

than a light-hour. Only a small mass of Fe is then required, which could have been 
readily produced by the star itself. 

o 

Subject headings: Gamma-rays: Bursts - X-rays - Cosmology: Miscellaneous 

6 

^ 1. Introduction 



Recently, Fe K-a and K-edge X-ray features have been reported with high significance in 
the afterglow of the bright burst GRB 991216 and GRB 000214 after about 1 day (Piro, et al, 
2000, Antonelli et al, 2000). A straightforward interpretation of this observation would imply a 
mass <; 10~ 2 — 1M of Fe at a distance of about one light-day, possibly due to a remnant of an 
explosive event or supernova which occurred days or weeks prior to the gamma-ray burst itself. 

Here, we suggest an alternative and perhaps less restrictive scenario: an extended, possibly 
magnetically dominated wind from a GRB impacting the expanding envelope of a massive 
progenitor star. This could be due either to a spinning-down millisecond super-pulsar or to a 
highly-magnetised torus around a black hole (e.g. Wheeler, et a/.2000), which could produce a 
luminosity that was still, one day after the original explosion, as high as L ~ 10 47 erg s -1 . This 
luminosity may not dominate the continuum afterglow; but we argue that it could be efficiently 
reprocessed, by a modest amount Mp e ~ 10 _8 Mq of material at distances ^ 10 13 cm, into an 
Fe line luminosity comparable to the observed value, together with a contribution to the X-ray 
continuum. Under this interpretation, the dominant continuum flux in the afterglow, even in the 
X-ray band, is still attributable to a standard decelerating blast wave. 



- 2 - 



2. Energy Input and Fe Line Strength 

The typical GRB model assumes that the energy input episode is brief, typically <5 1 — 10 2 
s, its energy and mass outflow being either a delta or a top-hat function. However, peculiarities 
in the early stages of some afterglows, e.g GRB 970508, have served as motivation for considering 
a more extended input period in which the energy or mass outflow rate may vary in time and 
the late energy input could exceed the prompt contribution (e.g. Rees & Meszaros 1998, Sari 
&: Meszaros 2000). Here we consider a related model, in which the power output continues at 
a diminishing rate, for a longer time of hours to days. The prolonged activity could arise if 
orbiting debris around a newly-formed black hole takes a long time to be completely swallowed, 
or if the central object becomes a fast-spinning neutron star, rather than a black hole. Collapsar 
or hypernova models (e.g. MacFadyen & Woosley, 1999; Paczyhski , 1998; Woosley 1993) or 
magnetar-like GRB models (Usov, 1994, Thompson, 1994, Wheeler et al. 2000) provide a natural 
scenario for a sudden burst followed by a more slowly decaying energy input. 

The power output would be primarily in a magnetically-driven relativistic wind (which, even 
during the later phases discussed here, would be hugely super-Eddington) . The generic magnetized 
outflow from a spinning compact object is 

L m ~ 2vrr 2 C (i? 2 /4vr) ~ 1.5 x 10 52 r 2 6 £ 2 5 erg s" 1 ~ 1.5 x 10 47 r 2 6 £ 2 2 . 5 erg s" 1 , (1) 

where the suffix denotes the B-field in gauss. The characteristic radius r Q = 10 6 r O 6 cm could be 
the inner radius of the accretion torus in a BH model, or, in a superpulsar model spinning with 
near-breakup angular velocity, the radius of the light cylinder. In a normal pulsar model, the 
field is assumed to maintain a steady value, and the luminosity declines as the spin rate slows 
down. However, during the early stages, B might decline more rapidly than the slowing-down 
timescale: the power output then declines in proportion to B 2 . For instance, if a torus is losing 
angular momentum due to magnetic torques, and is gradually draining into the hole, the magnetic 
field (and therefore the associated MHD torques and wind-driven energy losses) would gradually 
decline in step with the surviving mass. Alternatively, if the central object were a rapidly-spinning 
neutron star, the field may decline as uniform rotation is established and the higher moments 
decay away. 

There is no reason to expect that the field decay would follow a power law. However, just as 
an illustration, we note, from equation (1), that a luminosity decay L(t) cx i~ 125 , which would lead 
to 10 47 erg/s after one day, could be a consequence of a drop in B from from 10 15 G to 3 x 10 12 
G in a compact structure with stored energy of at least 10 52 ergs whose characteristic spin period 
remained constant (at a fraction of a millisecond) . 

A similar argument could be developed based on the concept of a-viscosity. For a hot dense 
torus around the BH resulting from collapse of the core of the progenitor star, the a-viscosity 
parameter is a ~ B 2 j '(4tt pv 2 ) ~ 0.05i? 2 5/ oj" 3 1 T 9 " 1 , and the viscous accretion time for a torus of 
outer radius 10 9 rg cm is t v i sc ~ (aQ)^ 1 ~ 10 r^ 2 B^ 2 5 pi^Tg s, while for a lower field of 3 x 10 12 G 
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it is tvisc ~ 1-5 rg B^2^ days for nT ~ constant, and the accretion of ^ 10 2 Mq in t ~ 1 day is 
sufficient to provide a characteristic L ~ 10 47 erg s _1 at a day. 

We envisage that the burst is triggered by the collapse of a core within a massive stellar 
envelope, as in the scenario developed by Woosley and his collaborators. A funnel along the 
rotation axis would have been blasted open during the 1-100 s duration of the original burst; 
it would subsequently enlarge, owing to the post-explosion expansion of the envelope of the 
progenitor star (e.g. Eichler & Levinson 1999; Woosley 1993). The ram pressure of the continuing 
MHD outflow would further enlarge the funnel and could, after one day, have expelled the envelope 
material from a region 10 13 cm across , even in the equatorial plane (this would require velocities 
of no more than 10 8 cm/s). 

The magnetised wind from the compact remnant (which we assume to be relativistic) would 
develop a stand-off shock before encountering the envelope material, and shocked relativistic 
plasma would be deflected along the funnel walls. In the absence of magnetic fields, the 
contact discontinuity between the shocked jet and stellar gas would have a tendency to develop 
Kelvin-Helmholz instabilities, which could lead to bulk heating of the stellar material. This 
would extend to large Thomson depths tt > 1, required to ensure a fluid behavior leading to 
non-magnetic K-H instabilities. However, magnetic fields are expected, since even if initially 
absent the instabilities would lead to mixing with the highly magnetized jet material. The K-H 
equivalent of highly relativistic MHD oblique shear instabilities is poorly understood, and it is 
unclear how deeply , if at all, such bulk heating would penetrate the thermal material outside the 
funnel wall. However, for energy deposition spread over a layer thick enough to have tt S> 1, the 
cooling rate (due to comptonisation, bremsstrahlung and recombination) would be high enough 
to reduce the temperature of the bulk-heated electrons and protons to the equivalent black body 
temperature, estimated as T e ~ T p ~ (Lm/aAirr^) 1 ^ 4 ~ fewlO 6 L^jr^^ 2 K. 

A more efficient heating mechanism of the stellar funnel wall gas is radiative heating, which 
would deposit energy within shallower layers with modest scattering optical depth. Non-thermal 
electrons are expected to be accelerated behind the standoff shock in the jet material; the 
transverse magnetic field strength (which decreases as 1/r in an outflowing wind) would be of 
order 10 4 G at 10 13 cm - strong enough to ensure that the shock-accelerated electrons cool 
promptly, yielding a power- law continuum extending into the X-ray band. Some of these X-rays 
would escape along the funnel, but at least half (the exact proportion depending on the geometry 
and flow pattern) would irradiate the material in the stellar envelope. For the high radiative 
efficiencies expected in relativistic shocks, the inward-directed radiative radiative heat flux would 
be comparable to that of a bulk heat flux from instabilities. However since the radiative flux is 
deposited in shallower layers containing less mass than the bulk heat flux, the radiatively heated 
shallow layer would be expected to be substantially hotter than a deeper bulk-heated region, its 
temperature being determined by photoionization equilibrium. 

The shocked relativistic material in the jet MHD flow exerts a pressure that is comparable 
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with the ram pressure of the pre-shocked outflow. This is aL m /47rr 2 c, where r = 10 13 ri3 is the 
distance of the funnel wall, and a 1 is a geometrical factor. This must be balanced by the 
thermal pressure 3n e kT in the photoionized layers of the stellar envelope forming the outside the 
wall of the funnel, giving a characteristic electron density 

n e = aL m /6Trr 2 ckT ~ 10 17 aL^r^Tf 1 cm -3 (2) 

The temperature parametrization used here is consistent with the range expected from 
photoionization equilibrium. The recombination time for hydrogenic Fe in the funnel walls 
photoionized by the non-thermal continuum produced by the stand-off shock is 

t rec = 6 x lO^Tg^n^ 1 = 6 x 10~ 6 a' 1 L^rl 3 T^ /2 s . (3) 

The ionization parameter is £ = (3L m /r 2 n e = lQ A (3a~ l T%. In this expression f3 < 1 is the ratio 
of ionizing to MHD luminosity: it is actually the product of two factors, namely the fraction of 
the MHD wind energy that is randomised in the stand-off shocks, and the fraction of that energy 
which goes into electrons that radiate in the X-ray band. For a large fraction of the Fe to be 
hydrogenic, £ must exceed 10 3 , and this condition would indeed be satisfied unless (3 were very 
small. The effective depth di to which Fe and other metals can be ionized is given by balancing 
the number rate of ionizing photons per square centimetre, /3L m /47rr 2 divided by their mean 
energy (about 10 kev), against the recombination rate nF e 5di/t rec . In this expression, <5 ^ 1 takes 
into account the recombination of other metals besides Fe. 

Standard calculations of photoionization of optically-thin slabs (e.g. Young, 1999) show that 
the equivalent width of the Fe K-alpha line, for solar abundances, is about 0.5 kev; the line is 
twice as strong, i.e. 1 kev equivalent width, if the Fe has ten times solar abundances (and even 
stronger for still higher enrichment). These results are applicable in the present context provided 
only that one further condition is satisfied: namely, that the ionizing photons encounter a Fe ion 
before being scattered by free electrons - provided, in other words, that tt = (TTdin e does not 
greatly exceed unity. We find that 

di ~ Zct rec x F \{(3/a5) ~ 6 x l0 9 ((3/a 2 5)T* /2 L^r 2 ^ 1 ^ cm (4) 

where xp e = 4 x 10^ 5 ^F e& is the Fe abundance normalized to the solar value. The Thomson 
depth through this layer tt = axdiUe ~ 4 x 10 2 {(3 / 'a5(Fe&)T^ 2 3 provided ((3 / a5CFeo) ^ 10~ 2 . 
The three parameters (3, a and 8 in this expression are all somewhat uncertain. However, the 
relevant material in the funnel and cavity walls comprises the innermost non-collapsed layers of 
the precursor star so it could well be greatly enriched in Fe (i.e. we might expect Cfsq to greatly 
exceed unity). 

Under the foregoing conditions the Fe K-q photon flux is about 0.1 of the X-ray continuum, 
and is 

N LFe ~ 10 54 L 47 /3 ph/s (5) 
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The signal observed t ~ 1.5 day after the GRB 991216 burst by Piro et al. (2000) corresponds, 
for an assumed distance D ~ 4 Gpc, to 6 x 10 52 ph/s, and that observed 0.5-1 day after GRB 
000214 by Antonelli et al (2000) at z ~ 0.47 corresponds to 3 x 10 51 ph/s. As is clear from the 
above discussion, a wind luminosity of 10 47 erg/s would be sufficient to yield the observed line 
signal provided that /?, the fraction of the power that goes into a photoionizing X-ray continuum, 
were not below 0.06. (3 is of course uncertain, but this number does not seem unreasonably high 
for the efficiency of particle acceleration by relativistic shocks, and suggests that our fiducial value 
of 10 erg s _1 for the overall luminosity need not be an overestimate. The continuum flux from 
GRB 991216 in the 1-10 kev band is observed to be 50-100 times stronger than the flux in the 
Fe line at 1.5 days. Since this factor would only be of order 10 for a photoionized slab with the 
properties envisaged here, this suggests that most of the continuum, even in the X-ray band, 
could still come from a standard afterglow model, with luminosity declining as a power-law in 
t, involving a decelerating blast wave. On the other hand, we cannot rule out the possibilility 
that much of the afterglow in all wavebands is due to a continuing power output from a compact 
remnant. In the latter case, the time-dependence could be more complicated (with possible rapid 
variability) and such effects should certainly be looked for. 

3. Discussion 

The total amount of Fe needed to explain the observed K-a line flux, arising in a thin layer 
of the funnel walls of a collapsar model, amounts to a very modest mass of Mp e ~ 1O~ 8 M , 
which could be Fe synthesized in the core. The Fe-enriched core material can easily reach a 
distance comparable to r ~ 10 13 cm in 1 day for an expansion velocity below the limit v ~ 10 9 
cm s _1 inferred by Piro et al. (2000) (c.f. also Antonelli et al 2000) from the line widths. 
Such subrelativistic velocities of the envelope material would arise naturally from low T shocks 
propagating through the star following the collapse, as well as from the burst explosion itself at 
larger angles from the burst jet axis. The natural progenitor scenario in which this can occur 
is a collapsar or hypernova model (Woosley, 1993, Paczynski, 1998) with "failed supernova" 
characteristics (i.e. little or no supernova display). 

The initial, energetic portion of the relativistic jet, with a typical burst duration of 1 — 10 
s, will rapidly expand beyond the stellar envelope, leading in the usual way to shocks and a 
decelerating blast wave. A continually decreasing fraction of energy, such as put out by a decaying 
magnetar, may continue being emitted for periods of a day or longer, and its reprocessing by the 
stellar envelope can be responsible for the observed Fe line emission in GRB 991216. Since the 
energy in this tail can decay faster than t -1 , the usual standard shock gamma-ray and afterglow 
scenario need not be affected, being determined by the first 1-10 s worth of the energy input. 

Finally, a few comments on how the present suggestion contrasts with other possible 
interpretation of X-ray lines in the afterglow. Because we invoke a continuing power output after 
the burst (in an MHD wind) the X-ray lines that are observed a day after the burst do not need to 
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come from material a light day (or more) in extent, as has been suggested earlier (Weth et al 2000, 
Bottcher, 2000). Such a large radius entails lower densities (and hence slower recombination rates), 
and thus requires a much larger total mass of Fe. But the most serious problem is that it may 
require a two-stage event: a supernova explosion, followed by a burst delayed by many days, so as 
to allow time for a sub-relativisic Fe-rich shell to have reached the requisite distance. Although 
such a model has been proposed (Stella & Vietri, 1998; Vietri, et al. 2000), detailed calculations 
suggest (cf Bottcher & Fryer 2000) that there are difficulties in understanding how an efficient 
burst could be generated afer such a long time-delay. Such models would have other potential 
observational consequences which may be used as a test, such as a reddening and flattening of the 
light curve at late times as inferred in some cases (e.g. Bloom et al., 1999, Galama et al, 1999). 

An alternative to the pre-ejection scenario has been proposed by Bottcher & Fryer (2000). 
They suggest that a very extended torus could be created during the inward-spiralling of a 
compact object through the envelope of a giant or supergiant - a process that could be a precursor 
of an exotic type of supernova. This supernova would not only generate the burst, but expel 
sub-relativistic material which is shock heated when it encounters to torus. The resultant thermal 
X-ray emission could display line features. Bottcher &: FryerUs suggestion is more attractive than 
the pre-explosion scenario, in that (as in the model proposed here) the X-ray lines come from a 
region much less than a light day across. However, it may not yield luminosities high enough to 
explain GRB 991216. More seriously, the shocked supernova ejecta which are postulated to emit 
the X-rays have a very large optical depth, which reduces the equivalent width of lines and would 
cause the ionization edge at 9.28kev to appear as an absorption rather than an emission feature. 
The large scattering optical depth may even (if it exceeds c/v) inhibit the escape of continuum 
radiation on a dynamical timescale. However, there could be particular geometries in which a 
modified scenario along these general lines could meet the constraints. 

Further data on X-ray spectral features from afterglows will surely offer important clues to 
the nature of the precursor star, and the compact object that triggers the burst. 
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